Time Reversal Breaking Superconducting State 
in the Phase Diagram of the Cuprates 
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We review and extend a previous study on the symmetry of the superconducting state, stimulated 
by recent tunneling and Andreev reflection measurements giving robust evidences for the existence 
of a d x 2_ y 2 + id xy order parameter in the overdoped regime of two different cuprates. Looking 
for a possible second-order phase transition from a standard d x 2_ y 2 to a mixed and time reversal 
breaking state, we confirm the results of our previous analysis on La2_ a; Sr a ;Cu04. In the case of 
Yi_j / Caj / Ba2Cu307_ a; as well, among all the allowed symmetries, the d x 2_ y 2 + id xy is the most 
favored one and the unconventional state is likely to occur in a small dome of the phase diagram 
located in the optimal-overdoped region and at very low temperatures. 
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I. INTRODUCTION 

Recent tunneling and Andreev reflection 
experiments[l, 2] on high-T c cuprate superconduc- 
tors (HTCS) have shown small deviations from the 
pure d x 2_ y 2 symmetry of the superconducting gap [3]. 
Deutscher and co-workers have determined the doping 
dependence of the zero-bias conductance peak in tunnel- 
ing experiments on (110) oriented Yi_j / Ca 2/ Ba2Cu307_ a; 
(YBCO) films and have shown that this peak, which is 
associated to the d x 2_ y 2 symmetry [4], is spontaneously 
splitted in some of the measured samples signaling the 
presence of a small secondary component of the bulk su- 
perconducting gap [5]. This secondary component, which 
removes the nodal lines of the d x 2_ y 2 primary one along 
the rX (diagonal) directions of the Brillouin zone, is 
compatible with a time-reversal breaking id xy gap func- 
tion. The splitting of the zero bias conductance peak is 
non-zero only for S > S op t (slightly overdoped samples), 
and vanishes for S < S op t (underdoped samples) . Heavily 
overdoped YBCO samples are not available. This sce- 
nario has been subsequently confirmed by Gonnclli and 
co-workers who performed a re-analysis of their Andreev 
reflection data on L^-xSr^CuC^ (LSCO)[2] originally 
interpreted with a d x 2_ y 2 + s symmetry. They showed 
that a better agreement with the experimental results is 



obtained assuming a d x 2_ 



id xy pairing symmetry[6] 



and found that the amplitude of the splitting is not 
monotonic with doping. They are able to perform 
measurements on heavily overdoped LSCO samples and 
find the sign for a closure of the d x 2_ y i + id xy region in 
the overdoped. 

In the past Balatsky[7] has studied the properties of 
the magnetic-field-induced transition to the d x 2_ y 2 +id xy 
state proposed by Laughlin[8] to explain the thermal con- 
ductivity plateaux in Bi2Sr2CaCu20s[9]. More recently 
a quantum critical point separating a d x 2_ y 2 + id xy from 
a pure d x 2_ y 2 phase has been suggested[10]. Recently 



we have given, in collaboration with other authors, a 
simple analysis of the gap equation based on general 
symmetry arguments that indicates the existence of a 
d x 2_ y 2 + id xy superconducting region in the phase dia- 
gram of the cupratesfll]. In accordance with the recent 
experimental results [1, 2] we have found that the uncon- 
ventional pairing state should be stable in the optimal- 
overdoped regime and at very low temperatures. 



II. SYMMETRY ANALYSIS AND BCS 
EQUATION 

In order to describe the transition between two 
different superconducting phases, as the one from a 
pure d x 2_ y 2 to a mixed d x 2_ y 2 + id xy state, we can 
use the classification of the superconducting classes 
based on the symmetry properties under crystal point 
transformations[12] . 

We start by imposing two constraints on the order pa- 
rameter: 

i. |A ff k| = |Ak| where gk denotes the vector obtained 

from k after the action of one of the Ng group 
operators O g . 

ii. Ak+K = Ak where K is a reciprocal lattice vector 
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The first condition assures that the dispersion = 
-y/£k + |A£| that enters the physical observables is a 
total symmetric function with respect to the lattice 
symmetry [13], while the second is simply to preserve the 
crystal periodicity. Given the strong evidence of singlet 
pairing in HTCS we require that 

in. A_k — Ak 

This further requirement imposes us to work with even 
functions of k, because of the antisymmetric spin part of 
the wave function. 

Due to the strong anisotropy of the copper oxides we 
take a two-dimensional square lattice described by the 
non abelian C^ v point group. This group has eight el- 
ements (Nq = 8): the identity operator E, one 2-fold 
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TABLE I: Civ characters table. In the last column on the 
right we have labeled the symmetry of the base functions for 
each representation. 



rotation C| and two 4-fold rotations C4 about the Z axis 
perpendicular to the XY lattice plane and four reflec- 
tions a v and a v > in the four vertical planes. 
We can write the order parameter as 



k 



(1) 



where the sum over \i runs over all the irreducible repre- 
sentations. Using standard methods of group theory[14] 
one can determine the expressions for the various 
in (1): 



(c) w M 



(k) 



(2) 



where the primed sum is restricted to non-equivalent sites 
of the lattice, that is sites r = (x, y) that cannot be ob- 
tained by means of the action of any of the O g operators 
on another site r'. We deal, in this way, only with gap 
amplitudes on a particular lattice distance while the sign 
relation between all the equivalent sites is entirely ac- 
counted for by the harmonics w['"' (k) 



Wr + (k) = cos k x x cos k y y + cos k x y cos k y x 



(k) = — sin k x x sin 



k y y + sin k x y sin k y x 

(3) 



< 

(k) = cos k x x cos k y y — cos k x y cos k y x 
(k) = — sin k x x sin k y y — sin k x y sin k y x 



For r = 0, the only non-zero harmonic belongs to the s + 
representation, and corresponds to the isotropic s-wave. 
For r = (±1, 0); (0, ±1) = 1 the s + and d + representa- 
tions are associated with 
waves, respectively. For r 



tions are associated with the extended s x 2 +y 2 and d x 2_ y 2 



vanishing harmonics are related to the s xy and d xy 



(±1,±1) = 2, the non- 
waves 

(s+ and d~ representations). In (3) we have not included 
the basis functions for the two-dimensional representa- 
tion E because requirement Hi. rules out odd functions 
of k. This is a peculiarity of the bi-dimensional square 
lattice. The T>^h group for instance, proper to a system 
in which the coupling between the different Cu-O planes 
is not negligible, possesses a two dimensional representa- 
tion which is even with respect to inversion and in that 
case it must be considered. 



It can be seen that, using the expansion defined by 
Eq.(l) and (2), the BCS equation becomes 



A<"> = -^(r^A^^^kV^VMk') 



where the kernel K (k) is defined as 



Jf(k) = 



tanh^ 




E.E^rW 
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(4) 



(5) 



and A r is a multiplicity factor equal to 1 for a generic 
r, 2 for a site on the axis or on the bisectors of the lat- 
tice and 8 for r = (0,0). It must be noted that the 
form (4) represents the mean-field equation for a sys- 
tem described by a generic static pairwise interaction 
between the electrons Vkk' (with V(r) the correspond- 
ing real space Fourier transform) and it is obtained with- 
out any previous decoupling of the interaction term in a 
given symmetry channel. It consists of an infinite dimen- 
sion homogeneous problem which, in general, will be of 
difficult solution. Nevertheless it can be observed that 
the appearance of a harmonic labeled by r is always con- 
trolled by the corresponding V(r). If the attractive po- 
tential is therefore a rapidly decreasing function of r we 
can restrict ourselves to the smallest distances bonds. We 
introduce the short-hand notations Vo, Vi, V2, . . . for the 
on-site, nearest- neighbor (nn), next-to-nearest-neighbor 
(nnn), ... components of V(r). Due to the repulsive 
character of the local interaction Vo in the HTCS, the 
isotropic s-wave is ruled out. We consider here attrac- 
tive (negative) nn (V\) and nnn (y 2 ) interactions. For the 
following discussion we do not take into account longer 
range terms. 

We have not yet made use explicitly of the first re- 
quirement on the invariance of the gap amplitude |Ak| 
under point group transformations. Its effect is to force 
the gap function to transform either like a single rep- 
resentation, or like a complex (time-reversal breaking) 

with A^ and A^ 



combination of the form A^ + iA^., 



belonging to two different representations of the point 
group. To describe this latter possibility, in the case of 
Vo > 0, V\ < 0, V2 < 0, we take the complex combina- 
tion A^w^(k) + iA^w^k), with both A^ and A^ real. 
Each gap parameter affects the other only through the 
quasi-particle spectrum and the BCS equations read 



.2 tanh(/?£k/2) 



2 tanh(/3£; k /2) 
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dk 
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(6) 
(7) 



where (3 is the inverse temperature and 
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with r\ = s + , d + , £ = s + , d~ (17 7^ C)- F° r typical values 
of the particle density in the HTCS, the Van Hove singu- 
larity (VHS) of the two-dimensional DOS strongly favors 
the d x 2_ y 2 symmetry, which thus represents the princi- 
pal component of the superconducting gap, in agreement 
with the phenomenology [15]. We take henceforth 77 = d + 
in Eqs. (6), (7). 

As we mentioned above the spectrum for the £ com- 
ponent is (pseudo)gapped by the principal d x 2_ y 2 gap. 
As a consequence a critical value | V£ r | for the onset of a 
mixed order parameter exists, as opposed to the case of 
the Cooper instability in the metallic phase. The critical 
coupling for both d xy and s xy secondary pairing is deter- 
mined by solving the BCS equations (6), (7) with = 
and £ = d~,s + respectively [11]. Numerical solutions 
show that the d x 2_ y 2 + id xy symmetry is much favored 
with respect to the d x 2_ y 2 + is xy as one expects from a 
qualitative analysis of Eq.(6),(7). Moreover, the contin- 
uous transition from d x 2_ y 2 to d x 2_ y 2 + is xy is generally 
prevented by a first-order transition from pure d x 2_ y 2 
to pure s xy , due to the strong competition between the 
two harmonics, which makes mixing unfavorable. This 
scenario does not find experimental support. The con- 
tinuous transition from d x 2_ y 2 to d x 2_ y 2 + id xy gap is 
therefore the most natural candidate if a pure d x 2_ y 2 
is to be modified by a small secondary component, as 
suggested by the experiments [1, 2]. As far as the depen- 
dence of the d xy critical coupling on the value of \V\\ is 
concerned it can be realized that the more | V\ | grows the 
smaller is the ratio |V^ cr |/| V\ | so as to favor the transition 
to the d x 2_ y 2 + id xy state in the intermediatc-to-strong 
coupling regime. 



III. RESULTS AND CONCLUSIONS 

In this section we show numerical solutions of 
Eq.(6),(7) obtained studying the dependence on hole con- 
centration (S) and on band parameters. In a previous 
work[ll] the phase diagram for the d x 2_ y 2 + id xy super- 
conducting region for the case of LSCO (closed Fermi 
surface) was derived. Here we want to show the case of 
bilayer compounds like Bi2Sr2CaCu208+5 (Bi2212) and 
YBCO (open Fermi surface). Our rigid band dispersion 
is thus £k = —2t(cosk x + cosk y ) + At' cos k x cos k y ± 
T;t±(cosk x — cosfcj,) 2 — /1 where t and t' are nn and 
nnn hoppings and fj, is the chemical potential. We take 
t = 200 meV, t' = 50 meV and t ± = 35 meV. For the 
case of Bi2212, ARPES measurements show a value for 
2f_L[16] that varies from 88 meV to 140 meV. Unfortu- 
nately we are not aware of similar data on YBCO, but the 
bare values of t± in Bi2212 and YBCO are similar, so that 
we can expect similar renormalized values. In addition 
we introduce a cut-off oj$ in the integrals which measures 
the typical energy scale of the attractive potential[17] . As 
shown in Fig. 1, for each doping 5, we fix V\ such that 
the critical temperature for d x 2_ y 2 pairing coincides with 
the pseudogap temperature T* (which we interpret as the 



temperature for Cooper pair formation without phase co- 
herence in the underdoped regime), and smoothly con- 
nects to T c around optimal doping. Then we fix a value 
for the ratio V2/V1 that never exceeds the physical thresh- 
old of 1. The connection between this value and the pres- 
ence of a quantum critical point located near optimum 
doping is given a previous work[ll]. 
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FIG. 1: hase diagram for d x 2_ y 2 and d x 2_ y 2 + id xy pairing in 
the HTCS. The parameters are appropriate for YBCO. 



We find that the mixed d x 2_ y 2 + id xy order parameter 
is then stabilized in a small dome slightly shifted in the 
overdoped region, 0.17 < S < 0.22. The d x 2_ y 2 + id xy 
transition temperature is always only a small fraction of 
the d x 2_ y 2 one and the same ratio is obtained between 
the gap amplitudes. In the extreme overdoped region 
the system approaches weak coupling determining the 
closure of the dome. On the other hand the vanishing of 
the d xy component in the underdoped region is due to 
band structure effects. 

We finally summarize the results of this analysis: 

1. As found in our previous work [11] a mixed order pa- 
rameter of the form d x 2_ y 2 + id xy can be stabilized in a 
small dome contained in the larger d x 2_ y 2 region. The 
d xy component does not grow indefinitely with increas- 
ing doping and it is always smaller than the d x 2_ y 2. A 
recent re-analysis of the LSCO data with a d x 2_ y 2 +id xy 
gap parameter [6], clearly displays the closure of the dome 
in the overdoped. The confirmation of our prediction in 
YBCO calls instead for further measurements at larger 
doping in the overdoped region. 

2. The intermediate-to-strong coupling regime is a nec- 
essary condition for the existence of the mixed order pa- 
rameter. 

3. d x 2_ y 2 + id xy pairing is more likely to occur in closed 
Fermi surface materials such as LSCO or electron doped 
materials, and in compounds where an interlayer hopping 
splits the VHS. 
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